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ABSTRACT 


Efforts  were  made  to  derive  the  design  criteria  of  surface  effect  vehicles  operated 
on  arctic  sea  !<_?.  Statistical  theories  were  developed  to  describe  trafficability  of  the 
vehicles  and  topography  of  the  sea  ice.  By  the  use  of  actual  sea  ice  surface  profiles 
obtained  by  an  aerial  laser  profiler,  the  usefulness  of  the  present  statistical  method 
was  demonstrated. 


MODEL  ANALYSIS  OF  VEHICLE  TRAFFICABILITY  WITH  APPLICATION 
TO  SURFACE  EFFECT  VEHICLES  ON  SEA  ICE  FIELDS 

Martin  Smith  and  Yoshisuke  Nakano 


I.  INTRODUCTION 

The  utility  of  a  surface  vehicle  is,  in  general,  a  function  of  the  vehicle’s  capabilities,  the 
terrain  in  which  it  is  to  be  operated,  the  operator's  skill,  and  the  navigational  aids  available. 
Further,  as  a  few  minutes  of  driving  in  a  snowstorm  makes  intuitively  clear,  even  the  most  detailed 
knowledge  of  the  above  variables  does  not  permit  a  sure  (non-statistical)  prediction  of  successful 
passage. 

A  straightforward,  effective  way  of  determining  trafficability  is  repeated  operation  in  repre¬ 
sentative  terrain.  If,  however,  one  is  interested  in  developing  a  vehicle  for  a  given  terrain,  experi¬ 
mental  studies  of  this  kind  can  be  an  expensive  and  time-consuming  method  of  optimizing  all  the 
test  vehicle’s  parameters. 

A  more  fruitful  technique,  particularly  in  the  early  design  stages,  is  to  develop  statistical 
models  which  enable  us  to  relate  various  parameters  of  vehicle  behavior  to  trafficability  in  a  given 
suite  of  terrains.  Such  models  may  enable  us  to  relate,  for  instance,  SEV  skirt  height  to  traffic- 
ability  in  typical  sea  ice  terrains.  The  extraction  of  statistical  predictions  of  vehicle  operation 
from  both  vehicle  parameters  and  terrain  data  is  the  object  of  the  present  study. 

Section  II,  together  with  the  Appendix,  develops  a  statistical  theory  of  trafficability  by  redu¬ 
cing  operator  influence  to  a  set  of  simple  decision  rules.  The  elimination  of  this  imponderable 
aliows  us  to  develop  a  concise  and  straightforward  theory  and  ensures  that  our  predictions  will  be 
lower  bounds  for  a  vehicle’s  behavior. 

Section  III  deals  with  the  extraction  of  pertinent  statistical  parameters  from  profilometer  data 
on  the  topography  of  sea  ice.  A  technique  based  upon  the  theory  of  random  functions  is  developed 
and  applied. 

Section  IV  combines  the  results  of  Sections  II  and  III  to  predict  a  relation  between  SEV  slope¬ 
climbing  capability  and  trafficability  in  sea  ice  terrains. 


II.  A  STATISTICAL  THEORY  OF  TRAFFICABILITY 
Statistical  descriptions  of  trafficability 

A  necessary  prelude  to  mathematical  models  of  trafficability  is  a  quantitative  notion  of  what 
trafficability  is.  What  statistical  parameters  of  a  vehicle’s  performance  are  of  interest  to  ua  in 
judging  its  utility  in  a  given  terrain?  A  number  of  candidates,  such  as  mean  velocity,  time  between 
breakdowns,  etc.,  immediately  present  themselves. 
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We  shall  consider,  here,  the  time  it  takes  a  vehicle  to  travel  a  given  distance  in  a  given  direc- 
tion.  From  this  definition  of  performance  the  following  statistics  naturally  follow. 

The  mean  time  of  arrival  is  the  time  we  would  "expect”  a  vehicle  to  require  to  make  the  jour¬ 
ney.  It  is  a  measure  of  the  vehicle's  operating  efficiency  in  a  given  terrain.  If  there  is  a  non-zero 
probability  of  a  vehicle’s  not  arriving  at  all  we  must  be  careful  in  our  definition  of  this  quantity. 

The  standard  deviation  of  the  time  of  arrival  is  another  quantity  of  interest.  It  gives  us  a 
measure  of  the  "spread"  of  arrival  times  we  should  anticipate  and  allows  us  to  set  "cutoff"  times 
after  which,  for  example,  search  and  rescue  operations  should  be  initiated. 

The  lateral  offset  is  a  measure  of  the  offset  of  the  vehicle  from  a  direct  path  from  origin  to 
destination.  Its  size  is  a  measure  of  the  vehicle’s  inability  to  pursue  a  direct  path  between  two 
points,  save  at  the  expense  of  increasing  travel  time. 

The  fatality  fate  is  the  probability  that  the  vehicle  will  not  successfully  complete  the  assigned 
mission.  It  is  a  measure  of  the  cost  of  a  mission  and  aids  the  rational  planning  of  rescue  capa¬ 
bilities. 


Suppose  that  we  conduct  a  set  of  N  experiments  in  which  vehicles  are  released  along  some 
given  initial  line  and  we  measure  the  time  t{  required  for  the  ith  vehicle  to  travel  a  distance  D 
away  from  the  line.  We  suppose,  further,  that  K  of  the  vehicles  never  arrive  so  that  tj  is  defined 
only  torN-K  values  of  i. 


The  mean  time  of  arrival,  then,  is  given  by 
1  * 

t  =  —  V  t, 

N-K* 


(1) 


where  1  denotes  summation  over  the  set  of  JV-K  successful  travelers, 
i 

The  standard  deviation  of  the  time  of  arrival  about  the  mean  is  given  by 


=r_L_  v  (t.  -  7)8> 

LN-K-l  Y  1  J 


(2) 


N-K-l  appears  in  the  denominator  instead  of  N-K  to  correct  for  our  use  of  t  estimated  from  the  same 
set  of  data. 

Let  v  be  the  mean  velocity  of  the  vehicles  given  by  Other  considerations.  Then 


t 


min 


D/v 


is  the  minimum  time  in  which  a  vehicle  could  make  the  trip,  and 


(3) 


is  the  "extra”  time  a  vehicle  spends  traveling.  If  we  assume  that  the  additional  time  is  spent  in 
traveling  in  a  random-walk  fashion  in  a  direction  normal  to  the  desired  course,  we  may  appeal  to 
random-walk  theory  to  estimate  how  far  to  the  side  of  a  straight-line  path  a  vehicle  is  likely  to 
emerge.  The  lateral  offset  LQ  is  given  by 


t0  -  v 


(5) 
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where  r  is  a  characteristic  length  for  the  time  interval  between  course  alterations,  r  may  be  thought 
of  as  the  typical  time  between  operator  decisions  regarding  the  vehicle's  course. 

The  fatality  rate  is  simply  given  by 

f  =  K/N  (6) 

and  is  the  fraction  of  vehicles  that  are  lost  in  transit. 


Statistical  modeling  of  trafflcabllity 

Bach  of  these  parameters  is  a  function  of  the  vehicle,  the  terrain,  the  operator,  and  the  infor¬ 
mation  available  to  the  operator  during  transit.  In  order  to  model  trafficability  usefully,  we  must 
reduce  this  dependence  to  a  small  number  of  relevant  parameters.  We  have  chosen  to  do  this  in 
the  manner  described  below. 


We  consider  the  vehicle's  universe  to  be  all  points  in  the  x-y  plane  given  by 


£  =  + nf 


(7) 


where  m  and  n  are  integers  and  ^and  fare  unit  vectors  in  the  x  and  y  directions.  Bach  such  point 
is,  at  any  given  time,  either  a  "passable”  point  or  an  "obstacle”  point.  We  may  define  a  function, 
t)  ( £ ),  over  the  plane  such  that 


) 


if  £  is  passable 
if  £  is  an  obstacle. 


(8) 


We  also  discretize  time  so  that  t  =  0,  1,  2 . *> . 

At  time  t  =  0,  the  vehicle  is  introduced  at  some  point  along  the  line  y  =  0.  The  vehicle  is 
assigned  the  mission  of  proceeding  "north"  (the  f  direction)  as  efficiently  as  possible  until  it 
achieves  the  line  y  =  IV,  at  which  time  it  ceases  movement.  During  a  unit  time  interval  the  vehicle 
may  move  one  point  east,  north,  or  west  of  its  previous  location. 

Let  £t  denote  a  vehicle's  position  at  time  t.  The  nature  of  the  vehicle's  next  mov^is  deter¬ 
mined  by  the  values  of  r\,  the  obstacle  function,  on  the  points  east,  north,  and  west  of  £t.  The 
precise  move  is  determined  by  application  of  the  following  decision  rules: 

1)  If  £t  +  f  is  a  passable  point,  the  vehicle  proceeds  to  it;  otherwise  it 

A  A 

2)  determines  if  both  of  £ t  t  Tare  passable  and,  if  so,  chooses  one  randomly  and  moves  to 
it;  otherwise  it 

a)  determines  if  one  of  £ ,  ±  f  is  passable  and,  if  so,  moves  to  it;  otherwise 

4)  the  vehicle  is  trapped  and  is  considered  a  fatality. 

Under  these  rules,  a  vehicle  is  just  barely  able  to  avoid  imminent  collision,  since  information 
is  available  only  about  the  immediately  adjacent  points.  Furthermore,  the  vehicle  cannot  reverse. 
Figure  1  portrays  the  decision  logic  graphically. 

A* 

We  see  that  if  r\  ( £)  =  0  for  all  £  the  vehicle  will  proceed  due  north  without  deflection  and 
will  achieye  its  mission  in  the  minimum  possible  time.  As  the  density  of  obstacle  points  increases 
from  zero,  the  vehicle’s  path  becomes  increasingly  circuitous  and  its  chance  of  entrapment  greater. 

We  have,  to  date,  considered  two  methods  for  assigning  values  to  the  obstacle  faction  i\.  The 
first  of  these  we  label  the  "scintillating”  or  "variable"  model.  At  each  time  t,  ( £ )  toe  a  given 
£  is  assigned  the  value  one  with  probability  p  and  the  value  zero  with  probability  1-p.  In  this 
model  we  have 
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West 


North 


East 


Decision  Logic 

o  Move  north  it  possible 

b.  Choose  omong  possible  eost  ond  west  moves  ot  rondom. 
c  Tropped 

Figure  l.  Driver’s  decision  logic. 


E  (l)  =  P  (9) 

where  E  is  the  mathematical  expectation  operator  over  either  space  or  time.  In  essence,  the  ob¬ 
stacle  space  is  recreated  at  each  time  t,  keeping  only  the  density  of  obstacles  constant.  We  have, 
to  date,  considered  only  spatially  and  temporally  homogeneous  distribution  of  obstacles. 

The  variable  model  simulates  the  effects  of  fog  patches,  drifting  snow,  etc.  Furthermore,  as 
discussed  in  the  Appendix  and  summarized  in  the  next  section,  we  have  been  able  to  extract  a 
number  of  useful  analytic  expressions  from  this  model. 

Because  the  behavior  of  the  obstacle  distribution  is  uncorrelated  in  time,  the  probability  that 
a  vehicle  will  be  able  to  make  any  particular  one  of  its  three  possible  moves,  or  none  of  them,  is 
independent  of  the  vehicle  s  history.  Consequently  each  move  of  a  vehicle  in  a  variable  apace  is 
one  sample  of  a  Markov  stochastic  process. 

The  second  method  we  have  used  for  assigning  7  differs  from  the  above  in  that  7  is  determined 
only  once,  at  time  t  =  0.  and  remains  fixed  thereafter.  We  have  labeled  this  the  ‘’fixed”  model. 

The  fixed  model  is  intended  to  simulate  the  effects  of  such  obstructions  as  terrain  features. 
Consequently,  we  consider  it  to  be  the  more  useful  of  the  two. 

The  passage  of  a  vehicle  through  any  given  fixed  space  is  also  a  Markov  process,  but  in  a 
less  useful  sense.  For  this  study  we  wish  to  find  the  statistics  associated  with  the  set  of  all 
fixed  grids  of  a  given  obstacle  density.  We  have  not  attempted  to  develop  any  analytic  results  for 
this  case  and  have  resorted  to  numerical  simulation. 

Figure  2  portrays  the  passage  of  three  vehicles  through  a  particular  fixed  space.  The  examples 
shown  are  illustrative  and  consist  of  a  vehicle  which  proceeds  directly  to  its  goal,  a  vehicle  which 
becomes  trapped,  and  a  vehicle  which  requires  three  extra  time  steps  to  complete  its  mission.  Note 
that  the  last  vehicle’s  second  move  could  have  been  into  a  cul-de-sac. 

Trafficability  in  a  variable  space  with  comparison  to  a  fixed  space 

We  present  here  some  analytic  and  numerical  results.  The  detailed  theory  of  trafficability  in 
a  variable  space  is  contained  in  the  Appendix. 

Let  p  denote  the  obstacle  density  and  K  tlie  number  of  north  moves  required  for  success.  We 
take  r  and  v  to  be  unity  ft  r  simplicity.  (This  scaling  can  always  be  achieved  by  an  appropriate 
choice  of  units  for  time  and  length.)  Then  for  a  variable  space  we  have 

K 

1  -  P  *  P3 
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O  =  \JK  (P-P3)  /  (1  -  p  +  p3) 

(10b) 

• 

Lq  =  2  v/K\  /  P~P —  ,  and 

V 1-P+P3 

(10c) 

0 

• 

>.[  1-P  )*. 

\ 1-P+P3 / 

(10d) 

0 

• 

We  may  also  define  the  dimensionless  incre- 

0 

mental  travel  time, 

0 

0 

t+  _  l~K 

K 

(11) 

x  which  is  the  fractional  increase  in  average  arri¬ 
val  time.  Clearly, 


Figure  2.  Typical  vehicle  paths  in 
fixed  obstacle  space. 


E  t+ 


7 +  _  P-P3 

1  -  p  +p3 


(12a) 


Since  t  is  independent  of  K ,  the  mean  relative  extra  travel  time  required  does  not  increase  with 
trip  length.  The  deviation  of  t+  about  its  mean  is  given  as 


E  [(t'-tVl 


-  i  h  M2' 

*  LPo  >o  /  J 


(12b) 


where  P(  are  given  by  eq  21-23  for  i  ■-=  0,  1,  2  in  the  Appendix  and  are  functions  of  p  only.  The 
deviation  of  t  decreases  as  K  increases,  since  each  vehicle’s  path  becomes  more  randomly  uniform 
with  increasing  trip  length. 

Table  I  shows  numerical  results  for  these  parameters  for  several  values  of  p.  The  columns 
labeled  PQ .  P,  and  P2  are  explained  in  the  Appendix.  K  has  been  taken  equal  to  24.  The  columns 
labeled  A.S.  represent  analytic  values;  those  labeled  M.C.  were  computed  by  Monte  Carlo  simula¬ 
tion.  We  have  shown  the  survival  rate  vSM,  defined  as 


SM  "  1_/  •  (13) 

The  rapid  decrease  of  SM  with  obstacle  density  suggests  that  most  results  of  practical  interest 
wiU  be  associated  with  small  values  of  p  ( <  0.1)  (or.  alternatively,  that  our  model  is  unrealistic). 

For  P  =  °-1*  we  can  8ee  that  the  8*tra  travel  time  (t  -K)  and  the  lateral  offset  are  not  large  com¬ 
pared  to  the  total  travel  time. 

We  define  the  penetration  depth  M*  as  the  distance  a  vehicle  has  a  70%  chance  of  ever  going 
in  a  variable  obstacle  field  of  a  given  density.  J If  is  given  by 

»*  -  log  (0.7) /log 

Table  II  shows  U  for  several  values  of  p.  Note  that  M*  decreases  very  rapidly  with  increasing  p. 
Table  III  permits  conversion  to  ovher  survival  rates.  Thus,  if  we  require  a  99%  chance  of  survival 

a  vehicle  may  proceed  only  0.0282  penetration  depths.  In  a  10%  obstacle  field,  then,  we  are  limited 
to  a  distance  of 
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Table  I.  Global  statistics  of  a  vehicle  in  a  variable  obstacle  space. 


Transition  parameter  SM 

P  p0  Pi  Pj  A.S.*  If. C .•* 


0.1 

.9989 

1.1086 

1.3636 

.9787 

.9739 

0.2 

.9901 

1.2264 

1.8229 

.7876 

.7870 

0.3 

.9629 

1.3244 

2.3894 

.4032 

.4036 

0.4 

.9036 

1.3609 

2.9667 

.0878 

.0877 

0.6 

.8000 

L2800 

3.3280 

.00472 

.00463 

•  A.S.  =  Analytical  Solution 
**  M.C.=  Monte  Carlo,  2.6  x  ID5  trials 


parameter  tor  K-24 


A.S. 

t 

If.C. 

A.S. 

a 

If.C. 

Lq 

26.637 

26.64 

1.720 

1.716 

3.848 

29.703 

29.70 

2.72^ 

2.664 

4.776 

33.012 

33.01 

3.761 

3.610 

6.004 

36.144 

36. 10 

4.909 

4.191 

6.970 

38.40 

38.66 

6.197 

4.846 

7.689 

Table  n.  Penetration  depth  is  a  variable 
obstacle  space. 


Table  III.  Snrvival  rate 
vs.  M/M*. 


P  IT* 


0.0 

oo 

0.1 

356. 

0.2 

36.5 

0.3 

9.46 

0.4 

3.53 

0.6 

LOO 

0.6 

0.826 

0.7 

0.467 

0.8 

0.281 

0.9 

0.189 

1.0 

0.000 

M/U'  SM 


.00281 

.999 

.0282 

.99 

.0864 

.97 

.144 

.95 

.296 

.90 

.626 

.80 

1.00 

.70 

1.94 

.50 

3.38 

.30 

4.61 

.20 

6.46 

.10 

8.40 

.05 

9.83 

.03 

12.91 

.01 

19.37 

.001 

Table  IV.  ComparlsoB  of  tine -of-ar rival  dlstribntloos 
for  variable  and  fixed  spaces.  If  -  24 


Time  ot 

arrival 

Variable 

p  =  0.02 
Fixed • 

R.E.7. 

Variable 

0.05 

Filed* 

R.E.7. 

Variable 

0.10 

Fixed* 

24 

.6157  ' 

.5892 

4.6 

.2903 

.2764 

6.0 

.0812 

.0786 

3.4 

25 

.2960 

.3169 

6.3 

.3469 

.3649 

2.6 

.1890 

.1941 

2.6 

26 

.0731 

.0754 

3.0 

.2203 

.2133 

3.3 

.23 24 

.2281 

1  q 

27 

.0131 

.0 136 

3.7 

.0965 

.0873 

9.4 

.2022 

.  1768 

14,5 

28 

.0017 

.0031 

44. 

.0329 

.0321 

2.6 

.1342 

.1094 

22.7 

29 

.0089 

.0126 

29. 

.0753 

.0626 

20.3 

30 

.0022 

.0065 

.0363 

.0366 

.9 

3 1 

.0164 

.0218 

29.6 

.0058 

.0134 

66.9 

*  Monte  Carlo,  2.6  x  105  trials  in  a  24  level  space. 
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Figure  3.  Survival  rate  vs.  trip  length 
for  various  obstacle  densities. 

0.0282  x  356=  10.04 

steps.  Figure  3  shows  the  rapid  decrease  in  survival  probability  SM  with  increasing  trip  length. 

Table  IV  shows  the  distribution  of  arrival  times  for  both  variable  and  fixed  spaces  of  depth 
24,  for  p  =  0.02,  0.05  and  0.10.  The  column  labeled  "R.E.%”  is  the  percent  relative  difference  Of 
the  two  adjacent  entries.  The  comparison  is,  we  think,  quite  satisfactory  for  those  time  intervals 
in  which  significant  arrival  probabilities  are  found. 

These  results  suggest  that  the  two  models  give  usefully  identical  results  for  low  obstacle 
densities.  Since  we  are,  in  practical  applications,  interested  only  in  low  densities,  we  can  apply 
the  above  analytic  results  to  either  model. 


III.  STATISTICAL  DESCRIPTION  OF  8EA  ICE 

Many  natural  phenomena  have  been  successfully  described  as  stochastic  processes  rather  than 
deterministic  ones.  One  of  the  earliest  and  most  typical  examples  of  such  a  situation  is  to  be 
found  in  the  theory  of  Brownian  motion,  where  each  coordinate  of  th  ;  Brownian  particle  is  a  random 
function  of  time. 

The  height  of  the  sea  ice  surface  can  also  be  considered  as  a  random  function  of  location  and 
time.  Hibler  and  LeSchack  (1970)  studied  sea  ice  surface  profiles  obtained  by  an  aerial  laser  pro¬ 
filer  approximately  200  miles  north  of  Point  Barrow  in  April  1968.  By  assuming  the  height  of  the 
sea  ice  surface  to  be  a  homogeneous  and  isotropic  random  function  of  location,  they  obtained  the 
correlation  function  and  power  spectral  density. 

We  extend  their  work  further  to  obtain  a  statistical  description  of  the  sea  ice  terrain,  which  can 
be  directly  connected  to  the  trafficability  of  the  SEV. 

Homogeneous,  isotropic,  and  normal  random  function 

We  introduce  a  random  function  H  (x)  defined  by 

H  (a)  =  height  of  sea  ice  surface 

where  x  =  a  point  (x<  x*2))  in  2-dimensional  Cartesian  coordinates. 
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In  addition  to  requiring  homogeneity  of  the  function  H(x),  we  introduce  two  assumptions;  namely, 
Assumption  1;  His  isotropic 
Assumption  2:  H  is  normally  distributed 

These  two  assumptions  have  not  been  validated;  however,  if  one  considers  a  large  area  o?  sea 
ice  terrain  these  assumptions  should  be  a  good  approximation.  We  define  the  mean  value  and 
correlation  functions  as 

£  H  (x)  =  m  (x)  mean  value  function  (15) 

£j[  H(x,)  -  MxjHfHfXj)  -  m  (*p  ]  1 

=  B  (Xj,  Xj)  correlation  function  (16) 

where  E  denotes  a  mathematic  expectation  operator. 

Since  H  (x)  is  homogeneous  and  isotropic,  B  is  a  function  only  of  the  distance  between  two 
points.  We  have,  then, 

®  (Xj ,  *j)  =  B  (r)  (17) 

where 


For  a  normal  random  function  H  (x),  the  mean  value  m  and  the  correlation  function  B  completely 
specify  H  (x),  i.e.  they  determine  all  the  multivariate  probability  functions  (Yaglom,  1962).  The 
multivariate  probability  density  functions  are  constructed  in  the  following  way  (Prabhu,  1965).  If 
the  covariance  matrix  T  with  element  =  B  (x, ,  x.)  is  known  and  is  nonsingular,  then  the  multi¬ 
variate  probability  density  function  /„  (h  r  hg . *n:  xr  x2 . xn)  is  given  as 

|y\|V4  (  i  n  i 

EZ  v*,-", '<*)-)>[  <«> 

where 

A  =  =  (Ajj)  , 

|A|  =  determinant  of  A, 
h  j  =  H(Xj),  and 
mi  =  E  H(Xj) . 

As  a  special  case,  the  two-point  probability  density  function  fg  (h hg ;  x,,  xg)  is  given  as 

|A|,/j  (  i  2  i 

exp  -  -  ^  £  *„  (»,  -  ■,)  (*j  -  »,) 

t  i.  j  =  1  ) 


(19) 
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where 


An, 

=  y 22/d 

*12 

II 

02 

*< 

II 

*22 

=  Yn/d 

d  = 

I'll  ^22  “ 

y  12 


Since  H  (x)  is  homogeneous  and  isotropic,  we  obtain 

i  r  i 


I g  (A  j*  ^2 '  O  - 


2 n  B(0)  \/£W 


exp 


2  B( 0)  (1  -  /92) 


{(hj-m)2  +  (h2  -  m)2-2/3(h1-m)(h2-m)j 


where 


/3  =  B(r)/B  (0)  . 


(20) 


Upper  bonnd  of  probability  for  finding  an  obstacle  derived 
from  two-point  probability  density  function 

Although  a  complete  quantitative  description  of  the  terrain  capabilities  of  an  SEV  has  not  yet 
been  obtained,  it  is  known  that  the  slope  of  the  surface  and  sudden  changes  in  sea  ice  height  are 
two  of  the  most  important  controlling  factors.  We  characterize  this  situation  as  follows.  Suppose 
we  are  operating  the  vehicle  at  a  certain  location,  and  we  are  interested  in  the  probability  of  finding 
an  obstacle  within  a  certain  distance  r  of  the  present  location,  x  =  x  r  We  begin  by  computing  the 
probability  that  H  (x  j)  =  h  t  and  H  (xg  =  x  {  +  r)  >  h  j  +  £;  namely, 


P*  IH  (x j)  =  /> j ,  and  H  (xg)  >  h  t  +  £  I 

no 

“  f  fg  j*  &  ^2 


■-  C  exp  (  -  D  (1  -  /S2)  (h  j  -  ro)2)  ll  -  <f>  (s)] 


(21) 


C  =  1/(2  ffB(0)vT^n 
D  =  1/  (2  B  (0)  (1  -/32)) 

y/TD  l(h1-m)(l-/8)+  <1 


s  = 


1 


dt 


where 
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For  SEV  operation  the  difference  in  the  heights  of  two  neighboring  locations  is  important 
rather  than  the  absolute  height  of  the  surface,  so  we  compute  the  following  probability: 

Pg  I H  (x2)  >  H  (x  t)  +  £ !  r  =  |  x  t  -  x2  j  t 

=  P2*  (r.O 

'  1  '  jmw  f  exp{-2iT(0) 


i  .  .  * 

*  [vrroy  1  “ '  ® (*  *  C'J "  *i  <82> 

The  function  P2  (r,  <£)  is  the  probability  that  height  of  the  two  points  separated  by  a  distance 
r  differs  more  than  £.  , 

Finally  we  compute  the  probability  P*  that  the  condition  H  (x  t)  >  H  (x  t)  +  £  is  ever  satisfied 
within  a  specified  distance  R.  Pg  is  given  as 


P* 

r3 


(£ 


■■>-T  / 


Pp  (f.  0  d  r 


(23) 


If  we  assume  that  the  SEV  is  designed  to  pass  over  an  ice  ridge  of  up  to  £(ft),  the  function 

S"  (C  R)  gives  an  upper  bound  on  the  probability  of  finding  an  obstacle  within  the  distance  R. 
w  (C  H)  can  also  be  interpreted  as  an  upper  bound  on  the  fraction  of  points  in  an  obstacle  state 
p  a?  described  in  Section  II,  if  the  distance  is  properly  chosen. 


IV.  DESIGN  CRITERIA  OF  A  SEV  DERIVED  FROM 
SEA  ICE  8URFACE  ROUGHNESS 

We  used  the  correlation  functions  obtained  by  Hibler  and  LeSchack  (1970)  from  natural  sea 
ice.  Figure  4  shows  three  samples  of  correlation  functions.  We  computed  the  values  of  P*  (£,  R) 
for  Sample  1  to  examine  the  general  behavior  of  Pg  .  The  results  of  this  computation  are  shown  in 
Figure  5.  in  which  the  values  of  P3*  are  plotted  as  a  function  of  R  with  £  being  a  parameter.  P  * 
naturally  increases  with  distance  and  decreases  rapidly  as  £  increases. 

As  we  mentioned  in  Section  III,  if  we  choose  the  distance  R  properly,  Pg  can  be  interpreted 
as  an  upper  bound  of  the  fraction  of  points  in  an  obstacle  state  p,  as  defined  in  Section  II.  We 
chose  the  distance  R  to  be  50  ft  as  a  trial  value  and  computed  Pg  (£,  50)  for  all  three  samples. 
The  results  of  this  computation  are  shown  in  Figure  6,  in  which  Pg  is  plotted  against  £.  It  is  easy 
to  see  from  this  figure  that  Sample  3  is  the  most  difficult  terrain  of  the  three.  We  note  that  even 
for  Sample  3,  if  ^ >_  2.0  (ft),  Pg  (£,  50)  is  less  than  10”^  or  1%. 

In  order  to  connect  the  results  of  Section  II  with  the  present  analysis,  we  plotted  expected 
relative  extra  travel  time  t*  versus  Pg  (£.50)  jn  Figure  7.  For  instance  if  the  total  distance  is 
1000  miles  and  Pg  =  0.1,  one  lias  to  travel  110  miles  extra.  .However,  if  the  SEV  can  pass  over  an 
ice  ridge  up  to  2.5  ft  in  height,  the  extra  distance  traveled  would  only  be  on  the  order  of  10  miles. 
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Figure  6.  P3*  (£  50)  vs.  £  tor  the  three  samples. 


Figure  7.  Expected  relative  extra  travel  time 
as  a  function  of  P*  ( £  R). 
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Figure  8.  Probability  of  being  trapped  every 
1000  miles  vs.  P*  (£  50). 


A  severe  restriction  is  imposed  by  the  fatality  rate.  In  Figure  8  we  plotted  the  values  of 
(1  -  &M),  which  is  probability  of  being  trapped,  for  every  1000  miles  versus  P  *  (/,  50).  where  S 
is  computed  as  follows  3  M 

103  miles 

SM  =  I  P3*  (C  50)  I 

.  p*  io5 
“  r3 


If  P3*  =  0.01, 
order  to  make  P* 


one  out  of  every  100  SEV’s  traveling  1000  miles  will  be  trapped  in  the  field. 
<  0.01  the  SEV  should  be  able  to  pass  over  about  2.5-ft  ice  ridges. 


In 


Although  our  present  model  needs  further  refinements  and  only  limited  information  on  natural 
sea  ice  is  available,  the  results  of  the  analysis  indicate  that  the  SEV  is  a  promising  operational 
vehicle  in  sea  ice  terrain,  in  view  of  the  Tact  that  a  small  prototype  SEV  is  able  to  pass  over  ridges 
up  to  3.0  ft  in  height  wit.ioit  difficulty.  Since  Held  studies  on  traffic  ability  are  often  expensive 
and  time-consuming,  this  model  analysis  may  provide  useful  guidance  to  designers. 
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APPENDIX:  STATISTICS  OP  A  VARIABLE  SPACE 

Lei  P(N )  denote  the  probability  that  a  vehicle  will  make  the  northward  transition  from  the  line 
y=i  to  line  y=/+ 1  after  exactly  N  moves.  We  assume  P(N)  is  independent  of  /,  time,  and  the 
vehicle's  east-west  position.  Clearly, 

P(N)  =  0  forN  <0  (1) 

We  define  the  moments  of  P(N),  Pq,  P j,  P2,  ....  by 

pL  -  t,  nm  »L  (2) 

N=  1 


Note  that  PQ  is  simply  the  probability  that  the  vehicle  will  ever  make  the  transition. 

Let  Q(,r)  (N)  denote  the  probability  that  a  vehicle  originating  at  y=0  at  time  zero  will  first 
achieve  the  line  y=K  at  exactly  time  K.  Clearly 

Q{IC)  (N)  =  0  for  N  <K  (3) 

We  extend  (3)  to  hold  for  negative  N. 

By  considering  all  the  ways  to  go  from  level  ff-1  to  level  K,  we  see  that  Q{K)  (N)  is  related 
to  (N)  by 

00 

o{K)m  =  £  q{K'1Hn-m)  p(U)  (4) 

M=  1 

P(M)  thus  plays  the  role  of  a  transfer  function. 

We  define  the  moments  of  (N)  as 


E 


N=  1 


Q(k)  (W)  flL 


(5) 


If  we  use  (4)  in  (5),  we  have 


00  00 


°L°  =  E  E  NL  QIk-»  (N -M)  P{M) 

N=  1  M=  1 


OO  OO 


E  E  W  -  M  +  *0L  15  (N  -  U)  P  (M) 

N= 1  M= 1 


E  E  11  (L\  (N  -  M)L~>  (N  -  M)  U1  P  (M) 

N=  1  M=  1  j=0  \;/ 


Preceding  page  blank 
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where^^is  the  binomial  coefficient.  We  may  reorder  the  sums  to  arrive  at 


°L°  =  E 


1=0 


CK 


(6) 


where 


ffLj  =  E  P  00  **j  ]C  1}  W  *  M)  (N  -  M)L~i 

M  1  N-  1 


or 


"lj  -  P)  oft11 


(7) 


Therefore 

0<L‘>  = 


(8) 


We  are  able,  therefore,  to  express  Dp  as  a  bilinear  combination  of  the  moments  of  Q^"1)  and 
P  of  order  less  than  or  equal  to  L.  For  example, 


Op  =  Po  °0*"l)  • 

Dp  =  PqD^'0  p  Pt  Dp_l>  , 

oP  =p0  °Tl)  +*8P1i>p«  ,  p2d^-d  , 

etc.  From  the  definition  it  is  clear  that 
Q(l)  (N)  =  P  (N) 

It  is  straightforward  to  show,  by  induction,  that  for  any  K  >  0  we  have 


(9) 

(10) 

(ID 


(12> 


op  =  (P0)K  . 

Dp  =  K  Pp1  P,,  and 


(13) 

(14) 


°P  =  *  ?$' “l  P2  +  If  df  -1)  P  f  Pj'2  . 


(15) 


The  mean  time  of  arrival  t  is  given  by 
f  =  Dp /D^. 


The  fatality  rate  f  is 
1  =  Dj*>  . 
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Tbs  standard  deviation  of  the  arrival  time  about  the  mean  is 

..kl  m»r 

Wv  J  ' 

The  lateral  offset  LQ  is 


From  (13)  -  (15)  we  have 


t  =  *  (P/P0) 

t  =  (P/ 


(16) 

(17) 

(18) 

(19) 


Our  remaining  task  is  to  relate  P(tf)  to  the  obstacle  density  p  for  a 
a  vehicle  has  just  achieved  y=/.  The  probability  that  it  will  be  able  to 
(1-P).  Therefore 


variable  space.  Suppose 
move  due  north  is  clearly 


P(D  =  (1-P)  . 


The  probability  that  it  will  not  be  able  to  move  north  but  will  be  able 
and  west  is  p  -  pd.  Then 


to  move  at  least 


one  of  east 


P(2)  =  (1-p)  (p  -p3) 

is  the  probability  of  moving  north  on  the  second  move.  In  general, 
P(N)  .  (1-p)  (p-p3)*"1,  N  >0  . 

We  insert  (20)  into  (2)  and  find. 


p  _  1~P 

1  -  p  +  p3  (21) 


(1-p  +  p3)8  (22) 

p  =  (1  -P)(l  +p  -p3) 

(1  -  p  +  p3)3  (221 

(21>  “  (23>  ^  be  substituted  directly  int0  (16>  -  (19)  to  yield  the  expression  given  in 


